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Abstract 

The versions of the self-consistent spin-wave theories (SSWT) of two- 
dimensional (2D) Heisenberg ferro- and antiferromagnets with a weak inter- 
layer coupling and/or magnetic anisotropy, that are based on the non- linear 
Dyson-Maleev, Schwinger, and combined boson-pseudofermion representa- 
tions, are analyzed. Analytical results for the temperature dependences of 
(sublattice) magnetization and short-range order parameter, and the critical 
points are obtained. The influence of external magnetic field is considered. 
Fluctuation corrections to SSWT are calculated within a random-phase ap- 
proximation which takes into account correctly leading and next-leading loga- 
rithmic singularities. These corrections are demonstrated to improve radically 
the agreement with experimental data on layered perovskites and other sys- 
tems. Thus an account of these fluctuations provides a quantitative theory of 
layered magnets. 
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I. INTRODUCTION 



Investigation of low-dimensional magnetism is an important branch of the modern solid 
state physics. Experimental interest in this problem is connected with the magnetic prop- 
erties of copper-oxide high- Tc superconductors, organic compounds, ferromagnetic films, 
multilayers and surfaces 

As stated by the Mermin- Wagner theorem, two-dimensional (2D) isotropic magnets pos- 
sess long-range order (LRO) only in the ground state. Unlike purely 2D Heisenberg magnets, 
real layered compounds have finite values of the magnetic ordering temperature Tm \J\ 
(J is the exchange integral) due to weak interlayer coupling and/or magnetic anisotropy. 
The smallness of transition temperature leads to some peculiar features of these systems. 
When crossing Tm, the short-range order (SRO) is not totally destroyed (in the 2D situa- 
tion it is maintained up to T ~ | J|), and a broad region above Tm with strong SRO exists. 
Corresponding experimental indications are provided by the data on elastic and inelastic 
neutron scattering: well-pronounced peaks of diffuse scattering were observed in La2Cu04 
0, Rb2MnF4 and K2NiF4 H, and well-defined spin waves in K2MnF4 up to T ~ 2Tn g. 

A great progress in the theory of the ground state and thermodynamic properties was 
made with the use of rigorous mathematical methods (quantum Monte-Carlo and renormal- 
ization group calculations). At the same time, simple approaches, which yield an analytical 
description of a wide range of physical properties, are very useful for practical purposes. At 
low temperatures (T -C Tm) the standard spin- wave theory works satisfactorily. At higher 
temperatures corrections owing to spin- wave interactions become important. These correc- 
tions were treated self-consistently many years ago for 3D Heisenberg model in Ref. P] The 
same results were obtained within a variational approach for isotropic and anisotropic 
Heisenberg magnets. 

For 2D magnets, close ideas were used recently by the "boson mean-field theory" 
which is based on the representation of spin operators through Schwinger bosons, and the 
"modified spin-wave theory" [|T^ based on the Dyson-Maleev (DM) representation. Note 
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that the former approach differs drastically from the standard mean-field approximation in 
the Heisenberg model: it takes into account spin-wave excitations and is highly sensitive to 
the space dimensionality. The results of these theories are in a good agreement with the 
scaling consideration |Tl| and experimental data on spin excitations in Cu02 -planes. Being 
generalized to quasi-2D case (see, for example, |T^-|T3[), these approaches lead to the same 



results as HH^]- The approaches of Refs. jSHlOl were also applied to frustrated 2D [|15Hig 
and 3D [^] antiferromagnets. 



In the approaches P,p!o| LRO is described in terms of boson condensation, see also Ref. 



112011 . To continue the region of applicability of the theory to disordered phase, a chemical 
potential of the Bose system (fictions magnetic field) is introduced at T > Tm, which is 
determined from the condition of vanishing of (staggered) magnetization. Introducing such 
a field can be more strictly justified within the projection operator technique . 



While the approach of Refs. corresponds to — oo limit of the generalized SU{N) 
Heisenberg model, the approach of Ref. |Ty] can be considered as the result of the self- 
consistent first-order l/S'-expansion, i.e. summation of all the bubble diagrams for the 
self-energy (see also Ref. 0). As argued in the present paper, this equivalence is preserved 
also between the first-order 1/N expansion and second-order self-consistent l/S'-expansion 
and seems to take place in all the orders of perturbation theories discussed. 

At the same time, above-discussed approaches (we refer them to as the self-consistent 
spin-wave theories, SSWT) turn out to have some shortcomings. The first one is mainly 
technical: the Bose condensation picture is inapplicable for anisotropic systems, since they 
have a gap in the excitation spectrum. Further, the SU{N) symmetry in this case is broken, 
so that the 1 /A^-expansion in the SU{N) model cannot be performed in principle. However, 
as it was mentioned, this expansion (and also description of LRO in terms of the Bose 
condensate) is not the only way to SSWT. 

The second shortcoming is much more essential. It concerns the description of ther- 
modynamics at temperatures that are comparable with Tm- In particular, the description 
of the (sublattice) magnetization curve near the ordering point is poor: the correspond- 



ing equation has two solutions, so that S{T) does not vanish continuously (see, e.g., Refs., 
PJI^). Besides that, the transition points are strongly overestimated. These drawbacks 
are due to that at sufficiently high temperatures the higher-order processes connected with 
dynamical interaction between spin waves, and also kinematical interaction should be taken 
into account. 

The kinematical interaction is important in a wide temperature region only for systems 
where Tm is not small in comparison with {JlS"^ (e.g., for 3D systems). This interaction 
can be, in principle, taken into account by the projection operator technique (see, e.g., Ref. 
pT|). However, this technique is rather complicated and is not convenient for practical 
purposes. Another way to obtain the corrections owing to the kinematical interaction is the 
use of the Baryakhtar-Krivoruchko- Jablonsky (BK J) representation [0,^ of spin operators 
via bosons and pseudofermions, which generalizes the DM representation (we do not know 
such a generalization for the Schwinger bosons). The introduction of pseudofermions gives in 
principle a possibility to exclude the contribution of unphysical states. These pseudofermions 
can be easily incorporated into the theory and, at least for 3D magnets, partially correct 
the above-mentioned drawbacks of the early versions of SSWT. 

For layered systems with Tm <^ | JjS*^ the kinematical interaction is less important (in 
fact, it works only in a narrow critical region near Tm), but higher-order (in the dynamic 
interaction) contributions should be included. As will be shown in the present paper, an 
infinite RPA-type series of diagrams are to be taken into account (as already mentioned, 
this is equivalent to the ffist-order 1/A^-expansion in the SU{N) model). Such a procedure 
permits to describe the "2D-like" Heisenberg regime [p4| , ^ where thermodynamics is de- 
termined by fluctuations of 2D Heisenberg nature. These results enable one to obtain the 
correct expression for Tm up to some non-singular constant in the denominator. At the same 
time, the true critical region, where the spin-wave picture of the spectrum is completely in- 
adequate, turns out to be very narrow in the layered systems. A satisfactory description of 
this region is possible within the 1/A^-expansion in the 0{N) model (see Refs. p^,|26[]). The 
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latter model is based on a fluctuation rather than spin-wave picture of excitation spectrum. 
This circumstance provides important advantages at high temperatures, but leads to some 
difficulties at the description of the low- and intermediate-temperature regions. Therefore 
the approach based on the SU{N) model (l/S'-expansion) is more appropriate at not too 
high temperatures. 

In the present paper we formulate a version of the SSWT, which is to a large measure free 
from above-mentioned shortcomings and is a good starting point for further improvements. 
To this end we (i) use the BKJ representation to obtain the correct description at not too 
low temperatures (ii) discard the condition S = and do not describe LRO in terms of 
Bose condensate, which permits to treat anisotropic systems. We also demonstrate (where 
possible) how our results can be obtained by the Schwinger boson method p. Further 
we calculate the corrections to SSWT for quasi-2D and easy-axis 2D magnets with small 
interlayer coupling or anisotropy using second-order spin-wave results in the self-consistent 
form. 

The plan of the present paper is as follows. In Sect.|I| we describe the representation 
of spin operators by Schwinger and Baryakhtar-Krivoruchko-Yablonsky ||22|. In Sects. ^ 
and we consider thermodynamics of quasi-2D and anisotropic 2D layered magnets within 
SSWT and construct an interpolation scheme between 2D and 3D cases. In Sect. 0we treat 
the problem of introducing magnetic field into SSWT and calculating magnetic susceptibility. 
In Sect. we investigate fluctuation corrections to the SSWT results, in particular to 
the (sublattice) magnetization and ordering temperature, and compare the results of our 
calculations with experimental data. 



II. BOSON REPRESENTATIONS OF THE SPIN OPERATORS 

We consider the anisotropic Heisenberg model 

H=-\i: j^A^^ -Ivi: J^,sts^ - D Y.{stf (1) 

ij ij i 
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where Jij are the exchange integrals, 77 > and D > are the two-site and single-site 
easy-axis anisotropy parameters. 



Consider first the Baryakhtar-Krivoruchko-Jablonsky representation [22,23 



St = V2Sb, , S^ = S- blk - (23 + l)c|c, (2) 



S- = V2S{bl - ^blblk) - ^-^^=^blch 



where 6|, bi are the Bose ideal magnon operators, and the auxiliary pseudofermion 

operators at the site i which take into account the kinematic interaction of spin waves. For 
the states \p) in the physical subspace (with the boson occupation numbers Ni < 2S and 
pseudofermion occupation numbers Fi = 0) we have c\p) = 0, and the representation (0) 
reduces to the standard DM representation. The states |mo) with A^, > 25, Fi = and 



with Fi = 1 are unphysical. As shown in Ref. ||22|, the partition function can be calculated 

as 



Z = Sp |exp ^—PH[b, b\ c, c^] — "ivr ^ c|c 



(3) 



where H[b, b\ c, c^] is the original spin Hamiltonian (|ip written through the Bose and Fermi 
operators according to (^. Analogous relations take place for the averages of spin operators. 
It can be proved that the contribution of states |mo) in (^) is exactly canceled by the 
contribution of the states 

Thus, unlike the DM representation, using the BKJ representation gives a possibility 
to exclude the contribution of the states with the boson occupation numbers Ni > 2S 
to thermodynamic quantities. It should be noted that this property relates to the exact 
Hamiltonian of boson-pseudofermion system H[b, b^ , c, c^] only and does not necessarily hold 
for its approximate expressions. However, one could expect that the introduction of the 
Fermi operators extends the region of applicability of approximate methods to not too low 
temperatures. 

The factor exp (^—in J2i c|ci) in (0) results in that the distribution function of the pseud- 
ofermions becomes —N{Ef) where N{E) = l/[exp{E/T) — 1] is the Bose function, Ef is the 



excitation energy for pseudofermions (as follows from the representation (0), the c-field has 
no dispersion). 

In the case of a two-sublattice antiferromagnet we separate the lattice into A and B 
sublattices. On the sublattice A we use the representation that is similar to 



S+ = y/2Sai , = S- aUi - {2S + l)clci, ieA (4) 



and on the sublattice B the "conjugate" representation: 



5+ = V2Sbl , = -S + h\h^ + {2S + l)4rf„ (5) 
S- = V2S{b, - ^b]bA) - '^^^^^^dtdA 



2S ' ' " V2S 

where a|, a^, and b], bi are the Bose operators, c|, q, and d\, di are the Fermi operators. 
Another useful representation of spin operators is the Schwinger-boson representation 



Si = XI 4a<^aa'Sia' (6) 



(TO-' 



where cr are the Pauli matrices, a, a' =1,1, so that 

St = ^(sltSiT - 4iSii), S+ = sl^Sii, S- = s^Si^ (7) 
The constraint condition 

4r^i^ + 4iSii = (8) 

should be satisfied at each lattice site. Since the phases of Sjf and Sii can be simultaneously 
changed, Sio- exp(i0j), this representation possesses a gauge symmetry. The Schwinger- 

boson representation can be simply related with the Holstein-Primakoff representation if we 
fix the gauge by the condition of hermiticity for one of the operators Sjo-, say, Sj|, i.e. = Sj|. 
Then we have from (M) 



Sit = ^25 - sJi^ii, (9) 



and substituting this into we obtain the Holstein-Primakoff representation. Thus the 
representations of the Schwinger bosons and by Holstein-Primakoff are equivalent. As well as 
for the BKJ representation, this equivalence can be violated in approximate treatments. Un- 
like the Holstein-Primakoff (or DM) representation, the Schwinger-boson representation can 
be easily generalized to an arbitrary number of boson "flavors" N > 2, and 1/A^-expansion 
can be developed. At the same time, there is no natural way to take into account "roughly" 
the kinematical interaction by introducing the Fermi operators into this representation . 

In the antiferromagnetic case we pass following to Ref. pl to the local coordinate system 
by the replacement 

at one of two sublattices. 



III. SSWT OF QUASI-2D MAGNETS 

A. Self-consistent approach within the BKJ representation 

In this section we consider the quasi-2D case with D = rjij = and the exchange integrals 
Jij = J for i,j being nearest neighbors in the same plane and Jij = J' for i,j in different 
planes. First we use the BKJ representation. In the ferromagnetic case the Heisenberg 
Hamiltonian (|l]) takes the form 

H = ~J2-J'^i [is - h\h - (23 + 1)c1q)(S - b]bj - (23 + l)c]c,) 

+25(61 - j^blblh)b, - 2(25 + l)blb,c^ (10) 
-^J2[blh + {23 + l)clci 

i 

To satisfy the condition 5* = in the paramagnetic phase we have introduced the Lagrange 
multiplier /i. This multiplier corresponds to the constraint of the total number of bosons and 
pseudofermions at T > Tc and plays a role of common chemical potential /i of the boson- 
pseudofermion system (for the pure boson system it was introduced in Refs. [^,0])- At 
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T < Tc we have fi = since no restriction of boson and pseudofermion occupation numbers 
is needed here. Introducing the chemical potential, which gives a possibility to continue the 
theory into the disordered phase, can be justified more strictly if one takes into account the 
kinematical interaction in a regular way |21]]. Since the magnon number is not conserved 
at T < Tc, the Bose condensation which takes place in Refs. P,pi]| does not occur in our 
approach. 



Further we perform decouplings of the quartic forms which occur in ([10|) . Introducing 
the averages 

j = S+{blh+sJ, i = S+{blb,+s,) (11) 
we derive the quadratic Hamiltonian of the mean-field approximation in the form 

H = J2 Js75 [bh - b\^,h + (25 + l)c\ci 

iS 

-^,Y.[blb^ + {2S+l)c^ (12) 

i 

where 75^^ = 7 and 75|| = 7'. From the definition of 7 (|ll|) one finds the system of self- 
consistent equations 

7 = + ^ A^k cos fca:, 'y' = S + ^Ni, COS k^, (13) 
k k 

which should be solved together with the condition 

S = S+{2S + l)N{Ej) - E (14) 

k 

where iVk = A^(-E'k) are the Bose occupation numbers, 

Ef = {2S + l){To-fi) (15) 
-Ek = To — Fk — /i 

are the pseudofermion excitation energy and the spin-wave spectrum respectively, 

Tk = 2 [I J|7(cos + cos ky) + \J'\j' cos k^] (16) 



Consider now the case of an antiferromagnet. Introducing the operators 

a,: i E A 

B. = { , a = { ^ (17) 




bl leB 
we derive 

Haf = I J|7 E [b\B. - BlsB, + (25 + (18) 
+ 1^17 E [b\B, - Bl,B, + {2S + l)ClC, 

i,<5|| 

-^^Y.[B\B, + {2S + l)ClC^^ 

i 

where 

^ = S + {aibi+sj_), y = S + {aibi+s^^) (19) 

Diagonahzing this Hamihonian one finds the self-consistent equations 

7 = ^ + E^cosA;,.coth^ (20) 

7' = S* + — — cos coth — ^ 
' u 2Ek 2T 



5 = (5 + 1/2) coth 4i - V ^ coth ^ 
^ ' ^ 2T Y 2:i 

where the antiferromagnetic spin-wave spectrum has the form 



^k = V(ro-/i)^-r2 (21) 

with Fk and Ef being the same as in the ferromagnetic case. 

For both ferro- and antiferromagnetic cases, the calculation of spin correlation functions 
PHI shows that is directly connected with the correlation length in the direction 6 by 
the relation 

^5' = 4^i^J\JM (22) 

The parameters 7 and 7' are also simply related to the spin correlation function at the 
nearest-neighbor sites, 
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\{SS^+s)\=7l (23) 
and therefore play a role of SRO parameters. For the total energy we readily obtain 

^ = -li: \j^.^s\i's = -mi' + (24) 

In the classical limit S ^ oo the SSWT equations are simplified. Supposing T ^ \J\S 
(Tm ~ \J\S^ in this case) the equations for both FM and AFM cases reduce to 

S/S = coth{Ef/2T) - I 5: ^ 

y^s + TY. ^ . (25) 

For T < Tm (/i = 0) the averaged (over nearest neighbors) SRO parameter 

7ef(T) = (4J7 + 2J'y)/^o (26) 
(but not the magnetization) satisfies the standard mean-field equation 

lei/S = (JoleiS/T) (27) 

where -800(2;) = cothx — 1/x is the classical Brillouin function (Langevin function). The 
temperature T* where 7ef(T*) = is higher than Tm, so that we have 7e/(7M) > and the 
behavior of 7e/ for T > Tj\/ is more complicated in comparison with (p^). 



B. Approximation of effective SRO parameter 

Equations (|13]), (p^ ) and (^) still demonstrate unphysical behavior of magnetization for 
T close to Tm at small enough J'/ J (see below). Introducing the pseudofermion field does 
not improve situation in this case: the transition temperature is already too small to be 
influenced by pseudofermion excitations with the energy of order of | J|. As discussed in the 
Introduction, the dynamical spin-wave interaction should be treated more correctly in such 
a situation. A rough solution of this problem can be achieved by the replacement 
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J2 Ji,i+slsiblbi - blbi+s) 7cf Ai+siblk - blbi+s) (28) 
5 s 

where 7ef is determined by (|26|). Then we obtain the spectrum 

= 7ef( Jo - Jq) " FM 



^q = V('^07ef - - ( Jq7ef)^ AFM (29) 

and the pseudofermion excitation energy 

= (25 + l)(7efJo-/^), (30) 

(here and hereafter we use the definition = J^slJsl exp(iq(5)). The SSWT equations take 
the form 

S = S+{2S + l)N{Ef)-Y.Ni, (31) 

k 

7ef = 5+^^JkiVk 

in the FM case and 

S=iS + 1/2) coth § - 7cf E ^ coth ^ (32) 

7ef = :S + ^V^coth^ 
Jo Y 2^k 2T 

in the AFM case. The approximation (^Sf ) is analogous to passing from the Hartree-Fock 
approximation to the local approximation in the spin-density functional method for the 
electron gas. As it is known, due to account of screening effects, such approximations can 
lead to improving the results and eliminating the unphysical peculiarities. Note that the 
same equations (|3TD and (^) were obtained earlier in Ref. However, when deriving 

these equations, the authors have used expressions for the spin Green's function which have 
incorrect u —>■ oo asymptotics. 

Another approach, which gives a possibility to improve the behavior of (sublattice) mag- 
netization near Tm, is based on a variational principle and is considered in Appendix A. 
It leads to the same spin- wave spectrum (p9|) , but the corresponding SSWT equations are 
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somewhat different from (^) and (|32|). However, numerical calculations shows that this 
difference is very small (several percents of magnetization value), and further we will refer 
to both these approaches as the approximation of effective SRO parameter. 



C. Temperature dependences of long- and short-range order parameters 

In the two-dimensional case J' = the spectrum E'k is independent of 7', and two 
remaining equations for S and 7 differ from those of approaches of Refs. P,[l2| only by 
the presence of pseudofermion distribution function N{Ef) which describes the kinematical 
interaction of spin waves. At T = we have fi = 0, and 7 = = S* in FM case and 
'J > S, S < S in AFM case, which corresponds to magnetic ordering in the ground state. 
As well as in Refs. equations ([T3|), (|14D and (^) do not have at T > solutions 



with yU = 0, S > since the integrals in equations (|T3|), (0) and ( pOf ) are logarithmically 
divergent in this case, and the only solution of these equations for J' = is S* = 0, yU < 0, 
which corresponds to a disordered phase. 

At low temperatures T ^ jJlS*^ we can neglect the pseudofermion contribution (i.e. 
kinematical interaction of spin waves) and we completely reproduce the results of Refs. 
fS HTOiP^ - In particular, the correlation length has the exponential dependence 



e = Cfexp(27rJ5VT) (FM) (33a) 
^ = Cfexp(2n\J\^oSo/T) (AFM) (33b) 



where 

So = S- 0.1971, 7o = ^ + 0.079 (34) 
are the 2D ground-state LRO and SRO parameters, Cj'^^ are the constants. The result 



was obtained earlier within the one- loop RG approach |rT|. With increasing T the 
role of kinematical interaction increases and for T ~ JS^ we cannot neglect N{Ef). The 
dependence ■j(T) for J' = is shown and compared with the result of approaches @JT^ in 
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Fig.|T] Unlike the approaches [|,|T3]; equations (0), (0) and ( pOD do not lead to the non- 



physical phase transition with vanishing of the SRO parameter, and the latter is finite at 
any temperatures. Note that for J' = the equations p9| ) and ( |32D give the same results. 

In the presence of interlayer coupling, the integrals in the SSWT equations (0), ( p^ and 
(PPP becomes convergent at finite T even at /i = 0. For not too high temperatures T < Tm 
(the ordering temperature Tm will be calculated below) these equations have the solution 
with S > 0, which corresponds to the ordered magnetic phase. For T > Tm we again have 
S = and /i < as well as in 2D case at finite T. 

Figs.|l]-§ show the results of the numerical solution of the equations of Sects. [Ill A| , [IIIB 



for different values of the interlayer coupling. In the three-dimensional case (J' = J) the 
(staggered) magnetization vanishes at Tc = 1.20 J (T/v = 1.33 |J|) that is approximately 
by 20% higher than the corresponding value obtained from the high-temperature series 
expansion. At the same time, the ratio T^/Tc = 1.20 is in agreement with the results of this 
expansion. The SRO parameter 7 demonstrates a sharp decrease in a narrow temperature 
region above Tm and then asymptotically goes to zero. One can see that 7afm > 7fm due to 
the quantum fluctuations. At the transition point we have 7^ = 7(T/v./) = 0.62 for FM case 
and 7c = 0.70 for AFM case. The dependences 'j{T)/S for 3D ferromagnets with different 
S are shown in Fig. |^ One can see that the value 7c/ S rapidly decreases with increasing S, 
reaching 7c = 0.39 at S* ^ 00. Thus at S = 1/2 strong quantum fluctuations are present 
even at T = Tm- 

Consider now the quasi-2D case < J'/ J < 1. At J'/ J < 0.4 equations ([T3|), (|14D and 
(pop still yield unphysical behavior of magnetization and SRO parameters for T close to Tm 
(as shown on Figs. |I] and for J'/ J = 0.3). At the same time, the approximation of single 
effective SRO parameter improves the behavior of magnetization for T close to Tm and 
provides a qualitatively correct description of thermodynamics at arbitrary temperatures. 
The price which we pay is overestimation of Tm even in comparison with the results of Eqs. 
(P!B|), (0) and (PPD, since the temperature dependence of the ratio of effective inter- and 
intralayer couplings (which is J'/ J for spectrum (^9]) ) is absent in the approximation used. 
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In particular, for J' / J — the results obtained within approximation ( pH]) are different 
from those of standard spin-wave theory only by quantum (ground-state) renormalization of 
7. Note that according to Fig.|I| with decreasing J' / J the size of the region with noticeable 
SRO increases. 

At small T — Tm we have — /i oc (T — TmY (see Fig. ^ for a ferromagnetic case, the same 
situation takes place in the AFM case) so that, according to (p2|), the critical exponent for 
the correlation length is z/ = 1. Since the magnetization changes linearly near Tm, we have 
also (3 = 1. The influence of higher-order terms in l/S* on these results is discussed in Sect. 
V|. Note that if we determine, following to Ref. [^, the critical exponent v from a not too 



narrow temperature interval near Tm , this becomes more close to the experimental value. 

At very low temperatures (T ^ |</'|'S') and arbitrary J'/ J the calculation can be per- 
formed analytically. The corrections to magnetization of a ferromagnet are proportional to 

2^3/2 

where C(3/2) is the Riemann zeta-function. At the same time, SRO parameters have a more 
weak T^/^-dependence 

Q / 7 \ 3/2 ^ s 5/2 



For J' = J this result corresponds to that of the Dyson theory ||28| to leading order in 1/S. 
For an antiferromagnet we have 

S = S, (38) 

24cJjJ'7o7o 



where 70, 7o and Sq are the zero-temperature values of corresponding parameters, c 



4J7o(2J7o + J'7o) is the spin- wave velocity. The corresponding temperature dependences 
of 7 and 7' are given by 
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In the case of small interlayer couplings J'/ J ^ 1 and higher temperatures, logarithmic 
singularities occur, and we can pick out them from the integrals in (p!3D, (jM]) and (^) in 



the same way as discussed in Ref. |jT3[ . In the quantum regime which takes place at not too 
low temperatures, where 

J'S<T<JS (FM) 
{JJ'f^S < T < \ J\S (AFM) (41) 



we obtain 



with 7 ~ 7o (the 2D values ( p4D can be used for 70 and Sq) and 7' being defined by the 
equation 

- T ( \ 

7' = 5o In — 1 (AFM). (43b) 

' " 47r|J|7 V 8JJ'77' / 

so that 7q = S'o- Note that in this case the infrared cutoff for the integrals over quasimomenta 
is 

{{T/JSf'^ (FM) 
90 = { (44) 
I T/c (AFM) 

(c = -\/8|</|7) rather than the boundary of the Brillouin zone. Since go 1? the continuum 
approximation for the excitation spectrum (and also interaction vertex) can be used in the 
quantum regime. Note that owing to the thermodynamic identity {dS/dT)s = {dS/dh)T 
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(with S being the entropy, h the magnetic field) the presence of T In T-terms in the magne- 
tization of a ferromagnet may be of interest in connection with the adiabatic coohng (see, 



e.g., Ref. H) 



For the critical temperatures in the regime (^T]) we obtain from (^21) the results 

47r|J|7,5o 

J-N 



HTysjj'^.Y,) 

with 7c = 7(Tj\/) ~ 7o and 7^ = '-^'{Tm) = 7j\,//47r| JI7. Comparing these results with the 
criteria of quantum regime (|4l| ) we obtain the condition of applicability of the results (^Sf) 
as 27r5' ^ ln(J/J'). It is important that 7^ <^ 7' and the interlayer coupling is strongly 
renormalized with the temperature. At the same time, only ground-state (quantum) renor- 
malizations are important for the intralayer coupling at \ J'\ ^ \ J\. 

In the case of large S (again supposing T ^ \J\S) we obtain for both ferro- and antifer- 
romagnet 

S = S--^\n^^ (46) 
An\J\S J'i ^ ' 

with 

7' = ^ - { In ^ - \\ (47) 

' 47r|J|5V J'i j ^ ' 

This leads to the expression for the critical temperature of a classical magnet with 1 -C 
ln( J/ J') < 27r ^ 

- ln(32J5/J'7^)- ^^^^ 
where 7^ = T/\,f/47r| JjS. As it should be, the critical temperature is the same for the classical 
ferro- and antiferromagnetic case. With the logarithmic accuracy we reproduce in this case 
the well-known results where iJS 1 (see, e.g., Ref. |^). Note that the factor of 32 which 
is often neglected leads to significant lowering of Tu as well as above- considered temperature 
dependence of 7'. 
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D. Mean-field Schwinger-boson approach. 



Similar results can be obtained within the Schwinger-boson representation. This is per- 
formed in the same way as in Refs. The Heisenberg Hamiltonian is written down 
in the form 



ij 



(49) 



where the chemical potential of bosons is introduced to take into account the constraint 
In the ferromagnetic case we subtract from the Hamiltonian (^Qf ) the term 



ij 



(50) 



to obtain 



<ij> i 



(51) 



where H = H — He, Tij = J2a ■sL-Sjo-, and : ... : stands for the normal ordering. Further the 
tilde at the Hamiltonian H will be dropped. Introducing the averages of the Bose operators 



7ii i'^ij) (^ij) 



(52) 



we derive the mean-field Hamiltonian 



Hmf = JiHij^ij - /^IK^It^^T + 4lSii] 

<ij> i 



(53) 



Such a procedure can be justified if we generalize the Schwinger-boson representation to the 
SU{N) model with arbitrary by introducing the operators (m = l...A^) and consider 
the limit ^ oo i]. 

In the quasi-2D case there are only two independent values of 'jij : 



7 i,j within the same plane 
7' otherwise 



(54) 



Introducing A = —fx — 7 Jq and passing to quasimomentum representation we obtain 

HmF = ^q*i<xSq<T (55) 

where i?q = A — Fq. Note that in the absence of external magnetic field the spectrum of 
bosons is doubly degenerate. The self-consistent equations have the form 



1 = Y1 COS k^, 7' = 51 COS h 

2S = J2N^^ (56) 

ko- 

As well as in Refs. [|],|12|, at low enough temperatures the Bose condensation takes place. 
Introducing external magnetic field (see Sect. |V|) removes the degeneracy of the boson 
spectrum, and only one of two bosons is condensed. Let A'kt (but not A^ki) contain the 
condensate contribution at /c — > 0: 

A^kT A^k + 2ns4o (57) 

where 2nB is the density of condensed bosons. Thus the self-consistent equations takes the 
same form as in the BKJ representation with S —>■ ub, N{Ef) = 0. 

In the antiferromagnetic case we subtract the term (pIlD from the Hamiltonian to obtain 
(cf. i) 

H = -^Yl Jii '■ -^Ij-^ij '■ -/^IZ(4tM + (58) 

<ij> i 

where Aij = sj^s'j^. Passing to the mean-field approximation we have 

Hmf = -^Y1 lijJiji-^ij + -Mj) - IK^T^iT + (59) 

<ij> i 

where 

= 7., = {A,) = {A\^) (60) 
Diagonalizing the Hamiltonian obtained one obtains 

Hmf = E ^q(«i«q + /?i/5q) (61) 
q 
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where = (A^ — r^)^/^. Thus the self-consistent equations take the form 



7 = E 7^ cos K{N^^ + iVki + 1), (62a) 



k 



r 



7' = E 77^ cos hiN^^ + iVki + 1) (62b) 



k 2^k 



r 



2^ = ETr(iVkT + iVki + l)-l (62c) 

k 

As well as in the ferromagnetic case, only A^kt contains the condensate contribution. Picking 
out this as 

A^kT/£^kT ^ NJE^ + riB (4o + 4q) (63) 

(Q = (vr, 71, 7i) is the wavevector of the antiferromagnetic structure) we get the SSWT 
equations (^2|) with coth(£'//T) = 1, S ^ ub- 

The corrections to above results can be obtained within the 1/A^-expansion in a general- 
ized Heisenberg SU{N) model (see, e.g., Refs. [p|,|32|-p5[] ) . As argued in the introduction (see 
also Sect. [V|) , the same results can be more easily obtained by higher-order l/S'-expansion. 
Thus the BKJ approach turns out to be more practical than the Schwinger-boson one. 



IV. SSWT OF THE EASY-AXIS 2D MAGNETS 

Consider now the 2D magnets with the easy-axis anisotropy. Besides the spin-wave exci- 
tations, the topological excitations (domain walls) contribute to thermodynamic quantities 
(see, e.g., discussion in Ref. |]3^). Such excitations cannot be taken into account in the 
approach under consideration. However, in the limit of small anisotropy 

D/|J|<l,r/<l (64) 

one can expect that the non-spin-wave excitations are important only in a narrow critical 
region. Outside this region thermodynamics can be described in terms of spin waves. Thus 
we restict ourselves to the case where (|6^) is satisfied. 
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Consider first the ferromagnetic case. Decoupling four-fold terms in the Hamiltonian we 
obtain 



(65) 



where 



^k = A-rk, Ef = {2S + l)\ 
X = Ml + vS)+D \{2S - 1) - A{blb,) 



(66) 



7 + V{bibi+s) 



and 



7 = ^ + {b%+s) 



(67) 



It should be noted that the expression for the excitation spectrum (^) is in fact the first- 
order 1/S expansion result. In particular, the spectrum (|66|) violates the requirement of 
vanishing of single-site anisotropy at S = 1/2 (this situation is discussed in Ref. 0). To 



correct this inconsistency we perform two replacements in the spectrum (|66D, which can be 
justified by calculating higher-order terms in 1/S: 



(25 - 1) - A{blk) -> (23 - 1) [l - 2{blb,)/s] ^ (23 - l)iS/Sf 

S - {blk+s) ^S\l- 2{blk)/S] \l + {b\h)/S - {b\h+s)/S\ ^ 



(68) 



Then the boson spectrum takes the form 



= 7(^0 - ^k) + JS^ - fi, 



(69) 



where 



A(r) = [{23 - 1)D/\JS\ + iJoS/J-f)v] (S/S)^ 



(70) 



is the dimensionless energy gap renormalized by spin-wave interactions. The system of the 
self-consistent equations reads 
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^ = S + ^J2JA (71) 
S = S -J2N^ + {2S + l)N{Ef) 

k 

In the antiferromagnetic case we obtain 

H = Y^ ^k(4«k + 4/5k) + Ef ^(4ck + 44) (72) 

k k 

with the spectrum 



Ek = VA2-r2, E/ = (25 + l)A (73) 
A = 7<^ + I^'S'IA - /i, Tk = 7Jk 

and A is the same as in ([70|). The system of the self-consistent equations takes the form 

7 = ^ + E ^ cos Z;;^. coth ^ (74) 

S-(S + i)coth|-i:^coth|i (75) 

Note that the proportionahty of the gap in the spin-wave spectrum to the squared sublattice 
magnetization was obtained earher within the renormahzed spin- wave theory |0,29|, which 



takes into account the influence of spin-wave interactions on the spectrum in a non-self- 



consistent way, and it is in agreement with the experimental data |30[ . 

Using the smallness of anisotropy and picking out the logarithmic singularities in the 
same way as in Sect. [lII Q we obtain 

Unlike the quasi-2D case, we have the unphysical result A{Tm) = because of the propor- 
tionality of the gap to {S/S)'^ (in fact a finite value of the gap at T = Tm should be caused 
by topological effects which are not taken into account). Thus we are unable to describe the 
dependence A(T) close to Tm- Denoting Ac = A(Tm) we have for the critical temperature 
at 2nS < ln(l/A) 
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T (77) 

ln(r/J5Ae)' ^ ' 

rj. 47r|J|5o7c 

J-N 



ln[TV8(J7c)2AJ 

In the case of large S we obtain for both ferro- and antiferromagnets 

— T 32 

This leads to the expression for the critical temperature of a classical magnet with 1 <^ 
ln(l/A) < 27rS' 

- ln(32/A,)- ^^^^ 

To leading logarithmic accuracy we can put Ac = A(0) in the above results. More correct 
calculation of Ac, as well as the corrections to the results ( [77D and (|79D will be obtained 
in Sect. 0. Note also that in the approximation A(T) = A(0), i.e. at neglecting the 
temperature dependence of the gap, we reproduce correctly the mean field result in the 
Ising limit, 

'S = SBsiJoSS/T) (80) 
where Bs{x) is the spin-S* Brillouin function. 

V. INFLUENCE OF THE EXTERNAL MAGNETIC FIELD AND THE 

MAGNETIC SUSCEPTIBILITY 

In this section we consider the influence of a weak external magnetic field h m a, ferro- 
magnet. This is described by the additional term in Hamiltonian, 

HH = -hY.St. (81) 

i 

The magnetic field results in an increase of magnetization, so that the total magnetization 
can be represented as 
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S — Ssp + Sind (82) 

where Ssp = S{h = 0) is the spontaneous magnetization, S'md is the field-induced part. 
Owing to the second term in the temperature dependence S(T) is changed: sharply 
decreasing in the vicinity of Tm, S nevertheless vanishes only in the limit T oo We 
consider a possible approach to the description of such a behavior in both versions of SSWT 
that are based on the Dyson-Maleev representation (or its generalization with the use of the 
BKJ representation) and Schwinger-boson representation. 

First we use the BKJ representation. The calculations, that are similar to described 
above, result in the equations (|13]) and (|1^ with the spectrum of spin waves 

Ek = l{Jo - Jk) + /i - /io (83) 

where /xq is the chemical potential in the absence of magnetic field: fiQ = at T < Tq and 
fj.o(T > Tc) is determined from the condition S'(T, /i = 0) = 0. 

Formally, the spectrum (^) has the same form as in the case of an anisotropic magnet 
(|55|) (we can associate with the anisotropy the effective "magnetic field" = JSA). How- 
ever, there is an important difference: in the case of the "true" magnetic field the chemical 
potential is taken at h = 0, so that the phase transition with vanishing S is absent (see 
below), while in the case of anisotropic magnet it should be determined in the presence 
of anisotropy field Ha, and S vanishes at Tc- However, at T ^ Tc this difference is not 
important (/i = /ig = in this region) and the magnetic anisotropy can be also described by 
introducing the temperature-dependent magnetic anisotropy field Ha- 

The temperature dependence of magnetization obtained by numerical solution of Eqs.(|13|) 
and (|1^) with the spectrum (^) is shown in Fig.^. At low temperatures T <^ Tc we have 
'S'sp ^ Sind and magnetization has mainly exchange origin. On the other hand, at T > Tc 
magnetization is entirely caused by infiuence of an external field and 

S^Xo'h,T:§>Tc (84) 

where 
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U^/h=o 4:T ^ smh^{Ej2T) T sinh^{Ef /2T) ^ ^ 

The first term in (pSj) differs from tlie result of tlie spin-wave tlieory by tlie form of tlie 
spectrum only, and the last term describes the correction owing to the kinematical inter- 
action. In a narrow region near Tq both contributions in (|82D are of the same order, and 
magnetization considerably differs from its zero-field value. 



It follows from (85) that Xo^ oc (T — Tm)~ so that the critical exponent is 7 = 2. Note 
that magnon-magnon interactions are taken into account in only by renormalization 
of single-particle spectrum. It is possible to improve result (^) by taking into account 
two-particle interactions in a RPA-type way, i.e. by considering the sum of the one-loop 
diagrams. This is performed in the next section. 

Now we consider the influence of external field in the Schwinger-boson representation. 
Carrying out the calculations similar to Sect. pi ^ we find for the boson spectrum 



Eka = la{Jo - Jk) - ^ha - fi (86) 

where = {slcr^icr)- The expression for the magnetization has the form 

5 = ^E(^kT-iVki)+ns (87) 

where we have taken into account the possibility of condensation of bosons with up "spins". 
There is also the condition of spin conservation at each site 

•^ = ^E(^kT + iVk|) + nB (88) 

At not too high temperatures T < Th, where Th is determined by the conditions 

^ = ^E(^kT + iVki), /i = -V2, (89) 
^ k 

the branch i?kt is gapless and > 0. AtT > both branches have a gap, and the condition 
(^) with ub = determines the common chemical potential. Thus the Schwinger-boson 
representation also allows us to describe the behavior of magnetization in the whole field 
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interval, and the expression (|87|) just describes magnetization as a sum of spontaneous and 
field-induced components. 

Up to now we have considered the small magnetic field values h J. In the opposite 
limit one can neglect the dispersion of boson spectrum (p3| ) and derive by using the BKJ 
representation the standard result 

S = SBs{Sh/T). (90) 



It should be noted that the correct result (|90[) is obtained only due to presence of pseud- 
ofermions, the Bose field alone leading to the unphysical phase transition with vanishing of 
magnetization at T ~ /i. 

A somewhat different situation takes place in the Schwinger boson representation. In 
the case h ^ J we have = and the equation for x = exp(— yu/T) has the form 

xcosh(V2T) - 1 



x2 + 1 -2xcosh(/i/2r) 
The solution to this equation reads 



S (91) 



^ = (1 + ^) cosh ^ + ^hs + If cosh^ A _ 45(5 + 1) (92) 



With the use of (0) we obtain the expression for the magnetization 

Thus in the limit of large magnetic fields the Schwinger-boson approach reproduces correct 
results only for 5 = 1/2. 



VI. FLUCTUATION CORRECTIONS TO SSWT FOR 2D AND QUASI-2D 

MAGNETS 

As already discussed, SSWT overestimates the value of Tm- In particular, for the simple 
cubic lattice the SSWT result for S* — > oo is Tm/S"^ = 1.803| J|. At the same time, the result 



of the spherical model (see, e.g., [^^ ) in this Imit reads 
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5*2 1 
3TAf Jq — J-k 



which coincides with the corresponding result of the Tyabhkov approximation |3^. One 



obtains from (Q) Tm/S^ = 1.319|J| which is close to the result of the high-temperature 
series expansion (see, e.g., Refs. ||49|j38|| ). As pointed in Sect. |Tll[ for S = 1/2 the value of 



Tm is overestimated by 1.2 times. 

In the quasi-2D case the formulas (^) and ( ^H]) (and the corresponding results of 2D 



case with small easy- axis anisotropy ( [77D and (|79D) coincide with the result of the Tyablikov 



approximation to logarithmic accuracy and thus seem to be correct. However, this accuracy 
is also insufiicient to treat experimental data (see detailed discussion in Ref. p4[) and the 
overestimation of Tm reaches 1.7 -r- 2.0 times for the quasi-2D case and nearly 1.5 times for 
anisotropic 2D case (the reason of weaker overestimation of Tm in the anisotropic case will 
be explained below). Thus in the quasi-2D magnets and 2D magnets with small easy axis 
anisotropy (in both cases Tm <^ l^l'S'^) the overestimation of Tm even higher then in 3D 
case. 

The values of critical exponents derived above {(3 = v = 1 and 7 = 2) are also in 
drastic discrepance with the molecular-field values {v = j3 = 1/2, 7 = 1), experimental 
data {u = 0.7, P = 0.33, 7 = 1.4) for isotropic magnets and exact values (z/ = 1, /? = 1/8, 
7 = 7/4) for easy-axis magnets, which are known from the Onsager solution of 2D Ising 
model. Thus SSWT describes poorly the critical behavior. 

At the same time, SSWT describes much better local properties (e.g., the pair spin 
correlation function at neighbor sites) than those determined by the scale of the correla- 
tion length. Indeed, at S = 1/2 the Tyablikov approximation yields the unphysical re- 
sult A£(Tc) = £{Tc) — £{S^) < [^. In the limit 5* ^ 00 this approximation gives 
A£^(Tc)/|£^(0)| = 0.6 which is also lower than the value which can be derived from the cal- 
culations in Sect.|ITl| (0.84). Besides that, the Tyablikov approximation implies a not quite 
correct form of the excitation spectrum at low temperatures. In particular, the spin-wave 
stiffness demonstrates the T^/^ dependence at low temperatures, instead of T^/^ one. 
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Generally speaking, the properties on the scales of order of correlation length cannot be 
treated correctly within one-particle picture, and the Tyablikov approximation gives only 
rough (but rather successful) description of these. A regular way of describing thermodynam- 
ics at not too low temperatures within spin-wave theory is to consider collective excitations 
rather than one-particle ones. For low-dimensional magnets with Tm ^ l^l'S'^, where large 



logarithms occur (see Sects, p^ll A| and [rV|) and fluctuations have 2D nature in a broad tem- 
perature region (except for the critical region), this can be performed analytically in a close 
analogy with the isotropic magnets of the dimensionality d = 2 + e (where j3 = 1 + 0{e), 



see, e.g., Ref. M). 

In this Section we take into account the interaction corrections to the SSWT results 
for the magnets with small inter layer coupling and/or anisotropy. Consider first the 2D 
Heisenberg magnet with the easy-axis anisotropy. In the ferromagnetic case we have 

H = Y.K^llK+\ ^(qi'q2;q3,q4)&i,&Ji2^qafeq4<5qi+q2,q3+q4 (95) 

q ^ qi...q4 

where 

E'^ = SiJo-J^) + \J\Sf 

<^(qi,q2;q3,q4) = Jcis + J^^ - Jm-^s - ^^-^i - -2|^l(qiq2 + /) (96) 

and 

f = i2S-l)D/\JS\ + iJo/J)vS (97) 

is the bare gap in the excitation spectrum. In the antiferromagnetic case, we use the oper- 
ators -Bq which are Fourier transformation of Bi of Eq. (|l7|) and satisfy 

aq= (5q + 5q+Q)/2 
bl^ = (fiq - i?q+Q)/2 (98) 

where Q = (vr, tt, tt) is the wavevector of the AFM structure. Then, up to some unimportant 
constant, we have the Hamiltonian of the same form (pSf), but for the operators i?q. Note 
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that in this case in (|95D has not the meaning of an excitation spectrum because of 



q 

non-Bose commutation relations for B^. 



[^q, K] = ^q,p + ^q.p+Q (99) 

The diagrams which give the first-order renormahzations of i?q and correspond to SSWT 
are shown in Fig. ||a (see, e.g., Ref. for the detailed description of this diagram tech- 
nique). Further on we suppose that all such renormahzations (which result in the replace- 
ments J — i> J'y/S and / — s> A in E^) are already performed and such diagrams can be 
omitted. 

To obtain the corrections to SSWT, higher-order diagrams should be considered. They 
lead to renormalization of one-particle energy (and occurrence of the damping) and also to 
vertex corrections. As discussed above, SSWT treats the excitation spectrum satisfactorily 
(this spectrum is already renormalized by first-order diagrams) . The calculations of damping 
of spin-waves, which occurs only in the second order of perturbation theory, shows that it 
is small in a broad temperature region Thus only vertex corrections should be taken 
into account. At not too low temperatures (T ^ | J|S'A) the RPA-type diagrams of Fig. ^ 
are most important since each loop contains a logarithmic divergence of the type ln(l/A). 
The integral equation for the vertex reads 

$(k, k - q; p - q, p) = if{k, k - q; p - q, p) (100) 

T ^ (^(k, k - q; s - q, s) 

2^ / 2 I A^^/ \2 I Ai ^(^' s - q; p - q, p) 



(^7)^r (^' + A)[(s-q)2 + A] 
(we have retained only the contribution of the modes with the Matsubara frequency u;„ = 0, 
which yields the logarithmic divergence, and dropped the terms with Un with simulta- 
neous cutting the summation over quasimomenta at the wavevector go which is determined 
by (|4^)). As can be seen from (^), the account of logarithmically divergent terms in the 
classical case can also be performed in the continuum approximation with = 32. The 
result of solution of the equation ( |100D (see Appendix B) in the 2D case reads 

2|J|7k(q-p) 



$(k, k - q; p - q, p) = = 



So - (T/27r I J|7) ln[go/ max( Ai/2 , g)] 
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-2|J|/ 



^ T q(q + k - p) q 



;ioi) 



7r|J|72 g2 ^1/2^ 

Note that the logarithmic corrections to the vertex in the isotropic case were obtained earher 
in Ref. [J^. For the static (staggered) non- uniform longitudinal susceptibility (for AFM case 



the shift q — s> q + Q is to be performed) we obtain from the diagrams of Fig.&; the result 



zz _ ~^qO 

A.Q 



l + {\J\^/2S)q^Xlf^ 



(T/27r|J|7So)ln[go/max(Ai/2,g)] 



(102) 



where 



T 

^90 - TJZyi X] 



(J7)2^(p' + A)[(p-q)2 + A] 
T/[27r(J7g)^]ln(gVA), » a ^^^^^ 
Xo = T/[47r(J7)'A], ^ A 

is the "bare" longitudinal susceptibility. Thus, as well as in RPA for itinerant magnets 



the spin susceptibility is enhanced by the interaction. It follows from the result ( |102|) that 
the excitation spectrum has different forms at small and large enough momenta: 

[ 25/(1 J|7g^), Wxl^-^S/S 

The first line corresponds to the standard spin-wave contribution (it is also subdivided in two 
cases as given by ( p.03| )). The second line corresponds to non-spin- wave regime: at ^ A 
one can neglect the anisotropy and xl^ oc is given, in particular, by the spherical model 
P3|| , which treats the spin excitations in essentially non-spin- wave way. Depending on the 
temperature value, three cases are possible. 

(i) low temperatures, T <^ Tm ~ 2tx\J\S'^ / \\i{ql/ /S). Then the second condition in ( |104[ ) 



cannot be satisfied and thus the excitations in the whole Brillouin zone have spin-wave 
nature. 
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(ii) intermediate temperatures, (5'/S')/ln(gQ/A) ^ T/27r| JjS*^ ^ S/S (T is of the same 
order as Tm)- Then at small enough q we still have Xq^ — Xqo^ but the second condition in 
( |104|) holds for large enough q where A exp(27r| J|7S'/T) <^ q^ < q^. 

(iii) critical region, T/2t{\J\S'^ ^ S/S (1 — T/Tm ^ !)• In this regime the first condition in 
( |104|) is satisfied only for q^ <^ A (hydrodynamic region) whereas at all other q the condition 
in the second line of Eq. ( |104[ ) is satisfied. 

The corrections to relative (sublattice) magnetization a = S/Sq (see diagrams of Fig.|^d) 
are given by 

^ V ^ I V $(k,k-q;k-q,k) 

\J\lSo ^k^ + A 2{\J\^rSo {k' + A)2[(k - q)2 + A] ^ ^ 



Integration leads to the result 

t 



In ^ + 4 In )— - 2(1 - Jf) + (t/a) 

Ao max((T, t) 



(106) 



where t = T/(27r| J|S'o7o)- The function $a takes into account the (unknown) non-singular 
contribution of non-RPA diagrams. Again, we have three temperature regions described 
above. In the region (i) only first term in the square brackets is to be taken into account 
and the magnetization demonstrates the spin-wave behavior ( [7S| ) and ( [75| ) for quantum 
and classical cases respectively. In the region (ii) all the terms, except for the last, are 
important, which leads to significant modification of the dependence S{T). The function $a 
in both regimes (i) and (ii) can be neglected and the result ( |106|) completely describes the 
behavior of magnetization in these two regimes. Finally, in the region (iii) the contribution 
of $a is of the same order as other terms in the square brackets. It should be noted that the 
factor of 4 before the second term in the square brackets is the sum of 2 which arises from the 
temperature renormalization of A^/^(T) oc AJ^^ max(a, t), and also a contribution of 2 arises 
from the vertex renormalization. Thus one can see that in the case of small anisotropy (the 
same situation takes place for small interlayer coupling, see below) the contribution from 
the renormalization of single-particle spectrum and interaction vertex are of the same order, 
so that SSWT is insufficient even outside the critical region. 

31 



For the Curie (Neel) temperatures we obtain the equations 

1 -1 



Tc = AnJS^ 
Tat = 47rJS'o7o 



\ Tc AtxJS^ ^ 
In ^„ . + 4 In h Cf 



JSA. 



c 



(FM) 



In 



T; 



N 



c^Ao 



+ 41n 



47r|J|^o7o 



+ C'af 



,32 ^, 47r|J|52 ^ ■ 
In — + 41n ' ' +Cei 

-^0 Af 



(AFM) 
(classical) 



(107a) 
(107b) 
(107c) 



with the constants Cf,af,c1 = — 2 — 4 ln2 + $^'^^''^^(00) which are still not determined within 
our approach. However, it is important that all the logarithmic terms are included in 
( |107|) and C give only a small contribution to above results. The gap at the ordering 
temperature, which remained indeterminate in Sect.|^, can be now estimated as A^ oc t^. 
The coefficient of proportionality is of order of unity and influences the constants C only. 

In the isotropic quasi-2D case the infrared cutoff for integrals over the Brillouin zone is 
J' / J rather than A. Then we obtain in the same way 



a 



max(cr, t) 



2(l-a) + <l>i, {t/a) 



(108) 



and 



1 -1 



Tat = An\ J|S'5'o7o 
Tm = ^ttJS^ 



J'\S 
In 



'c 



-'at 



8|JJ'|7o7o 
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In — + 3 In 

Aq Tm 



+ 3 In 



47r| 715070 



T 



N 



AF 



(FM) 

-1 



(AFM) 



(classical) 



(109a) 
(109b) 
(109c) 



In this case we have 7'(T) oc 7Qmax(a, t) which leads to that the coefficient at the second 
term in the square brackets is 3 (instead of 4 in the anisotropic case). This is why the 
interaction corrections are weaker in the anisotropic case: within SSWT the above-mentioned 
coefficient is 1 in the quasi-2D case (which is 3 times smaller than the correct value) and 2 
in the anisotropic case (only 2 times smaller than the correct value). Note that the results 
( |108| ), ( |109| ) are valid for all the four combinations of the signs of intra- and inter-plane 
exchange integrals (for mixed combinations, FM and AFM denote the type of the in-plane 
ordering) . 
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The same results ( p.06| )- ([T09|) were obtained within the RG approach in Ref. (Note 
that different sign at the third term of square brackets of ( |106[ ) and ( p.08| ) is the misprint 
of this paper). Derivation of general expressions for the case where both interlayer coupling 



and anisotropy are of the same order can be also found in Ref. p5|. 



With neglect of the functions $a(a;) and $ic(a;), Eqs.( |106| ) and ( |108| ) still yield unphysical 



behavior near Tm- The point T*, where the derivative dS/dT diverges, can be determined 
from the condition 

3/2 quasi-2D 

a(f)/r^<{ (110) 
2 easy- axis 2D 



which should be used together with (|106|) or (|108|) ; t* is the value of t corresponding to T*. 
The condition 

The functions and $ic(a;) describe the crossover from an isotropic 2D Heisenberg 

to 2D Ising and 3D Heisenberg behavior respectively. As discussed above, these functions 
give considerable contributions in the crossover region between regimes (ii) and (iii) and 
in the critical regime (iii), where essentially non-spin- wave excitations should be taken into 
account. An account of these functions results in slight decreasing the temperature T* in 
comparison with that given by (|110|) , and T* becomes the temperature of a rapid decrease 
of S (in fact, the characteristic temperature of a crossover). For a quantum antiferromagnet. 



the calculation of $ic(a;) can be performed within the 1/A^ expansion in 0{N) model pi 



For an arbitrary x = t/a, the result of this calculation is very cumbersome. In the critical 



region (x ^ 1) it provides the correct critical behavior p4 

^Neel 



a' 



47r|J|S'o7o 



r 1 / X n 2/33 



1 ^ 



1111 



-1 — ^0 V TnccI ^ 

with Ao ~ 0.9635 and (3^ = {1 - 8/7r^N)/2 ^ 0.36. The value of C'^p obtained by this 
expansion is very small, C^p = —0.0660. Other critical exponents can also be calculated 
within the 1/A^ expansion in 0{N) model (see, e.g., Ref. 



1/3 = 1- 32/3n^N ~ 0.64, 73 = 2(1 - U/tt^N) ~ 1.21 
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(note that the scahng relations are shghtly violated because of approximate character of this 
expansion for = 3). Thus the results of the spherical model for the critical exponents 
above Tm (see, e.g., Ref. |^) become radically improved. In particular, the fluctuations 
correct the critical behavior of magnetic susceptibility. 

For practical purposes, it is useful to have simple interpolation expressions for the func- 
tions $(x), which enable one to describe the crossover temperature region. Taking into 
account the closeness of T* to Tm which is given by ( |107|) and ( |109| ) and using ( |110| ), we 
can write down the simplest expressions for $(a;) in the form: 

$a'^"'^'(^) = ^|=(Cf,af,c1 - 2 + 81n2) 
va; + 1 

$r/"'^'(x) = -^=^(C^,,p^,i - 1 + 31n3) (112) 

(x < 1). The constants Cf,af,c1 and Cp^j can be in principle obtained from numerical cal- 
culations or by comparing with experimental data (see below). However, one should expect 
that they are small enough and can be neglected. 



VII. COMPARISON WITH EXPERIMENTAL DATA 

To discuss the experimental situation, we consider flrst the compounds with layered 
perovskite structure. The parameters used are given by Table 1. The experimental val- 
ues of transition temperatures are also given and compared with the theoretical one (for 
experimental data see Ref. |^9| and references therein, and Ref. for La2Cu04). 



Table 1. The experimental parameters and ordering temperatures of layered magnets 
and the corresponding calculated values of Tm in the standard spin- wave-theory (SWT), 
SSWT and RPA (in brackets - with account of the constant Caf = —0.7). 
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Compound 


S 


J,K 


J',K 


Ao 




T^SSWT TV- 
-'m 






La2Cu04 


1/2 


1600 


0.8 


^0 


672 


537 


343 


325 


K2NiF4 


1 


102 


^0 


0.0088 


160 


125 


90.0 (97.0) 


97.1 


Rb2NiF4 


1 


82 


^0 


0.046 


180 


118 


88.4 (95.0) 


94.5 


K2MnF4 


5/2 


8.4 


^0 


0.015 


74.8 


52.1 


42.7 (45.1) 


42.1 


CrBrs 


3/2 


12.38 


1.0 


0.024 


79.2 


51.2 


39.0 


40.0 



The values of J' and Aq = A(0) are obtained from the low-temperature behavior of the 
sublattice magnetization. This procedure gives a possibility to determine correctly the 
parameters since the results obtained in Sects. |III Q and work well at low temperatures 
(in particular they give correct results for ground-state renormalizations) . Note that for the 
anisotropic perovskites the parameters obtained are also in agreement with the experimental 
data on the spin- wave spectrum . It should be stressed that the experimentally observable 
gap in the spin-wave spectrum is A(T) whereas D,-)] plays the role of the bare parameters. 
The same situation takes place for the in the quasi-2D case where J'~f/S and J''y'/S are 
experimentally observable rather than the bare parameters J and J'. Since in the systems 
under consideration the parameters 7' and A have strong temperature dependence (see 
Sects. |111 C| and [TV] ) it is important to take into account this dependence when treating the 
experimental data. 

One can see from the Table 1 that for all the systems the estimated values of transition 
temperatures are close to the experimental results (for La2Cu04 the experimental data on 
J' are contradictionary; one of the possibilities to improve the agreement with experimental 
data is introducing small easy-axis anisotropy ^5]). At the same time, using SSWT without 
fluctuation corrections overestimates T/v by about 1.7 times, although improves somewhat 
the results of standard spin-wave theory. For the anisotropic compounds T^r is slightly 
underestimated. This may be due to two reasons: non-zero value of Caf in this case (Caf — 
—0.7 is obtained by best fit to the experimental data) and (less important) small interlayer 
coupling which also increases the transition temperature. 
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The temperature dependence of the sublattice magnetization for K2NiF4 is shown and 
compared with different theoretical results in Fig. 0. The regimes (i),(ii) correspond to 
T < 80K where the RPA (RG) result is in good agreement with the experimental data. 
It can be seen also that the account of the function given by ( |112| ) improves 

considerably the agreement in the crossover temperature region. The 1/A^-expansion in the 
0{N) model also gives satisfactory description of this region, but does not describe 
correctly low enough temperatures, since, as discussed in the Introduction, it implies an 
essentially non-spin-wave picture of the excitation spectrum. 

The parameter values and Curie temperature for the ferromagnetic compound CrBra 



(see Ref. for experimental data) are also presented in Table 1. Here both interplane 
coupling and anisotropy are important, and we obtain (see also Ref. p5| ) 



-M 



47r| J|S'o7o 



X <^ ln[2g^/(Ae + 2«, + JA^ + 4acA,)] + 2 In 



4vr|J|5o7o 



+ C 



-M 



where 



A 



a 



(113) 



a, = J'lUiJlc) = Tm\J'\/ [47r(J7o)' 

2 



Ae = Ac 



Tjv// (47175070) 



(114) 



The (unknown) function C{A/a) satisfies C(0) = C and C(oo) = C where C and C 
are defined above. Since the in-plane lattice structure is non-square, we have used the 
effective value of in-plane exchange integral determined in the continuum limit from the 
excitation spectrum as -Eq ^ 7</ef(9x + 9'v) + ^('^'5 J^)- One can see that the agreement with 



experimental Tc is excellent. 



VIII. CONCLUSIONS 



In the present paper we have investigated in detail the capabilities of the self-consistent 
spin-wave theory, which is based on boson representations for spin operators, for description 

36 



of layered magnets. To improve the SSWT at not too low temperatures, we have introduced 
a pseudofermion field. For magnets with low transition point, Tm ^ l^l'S'^, analytical 
results were obtained. These results have different forms in quantum (T ^ 1^15') and 
classical (T ^ l^l'S') regimes. In the quantum case the magnetization demonstrates in some 
temperature region the TlnT-behavior. The proposed version of SSWT gives a qualitative 
(and at not too high temperatures - quantitative) description of thermodynamics of layered 
magnets. An important advantage of SSWT (in comparison with the methods that are 
based on investigation of continuum models, e.g., non-linear sigma model) is a possibility to 
describe the short-range order above the transition point. 

At the same time, even in the case of layered magnets SSWT is unable to treat quantita- 
tively the transition points and thermodynamics at high enough temperatures (T > O.STm)- 
We have performed a systematic treatment of corrections to SSWT in the case Tm ^ I'^I'S'^, 
which is based on summation of higher order 1 / S'-terms. The inclusion of RPA corrections 
(which permits to take into account next-leading logarithmic singularities) yields an excel- 
lent description of the behavior of (sublattice) magnetization at arbitrary T < Tm except 
for a narrow critical region, where an account of non-spin-wave excitations is required. The 
approach used is somewhat reminiscent of the theory of itinerant magnets As well as 
in the latter case, the fluctuation corrections within the RPA approximation lead to sig- 
nificant lowering of the transition temperature and improve radically the agreement with 
experimental data. The simple analytical results obtained give a quantitative description 
of the magnetization behavior practically up to Tm- At the same time, the consideration 
of the critical region (e.g., correct calculation of critical exponents) requires an account 
of essentially non-spin-wave excitations. For quasi-2D isotropic magnets this can be per- 
formed within the 1/A^-expansion in 0{N) model ||2^. A description of the critical region 



for magnets with the easy-axis anisotropy, where the topological (domain-wall) excitations 
are present, is still an open problem. 

A regular calculation of higher-order corrections in the 3D case, where the kinematical 
spin- wave interaction should be also taken into account and the logarithmic terms are absent, 

37 



is also the matter for further investigations. The same is vahd for magnets with essential 
role of topological excitations (easy-plane systems, antiferromagnetic half-integer spin chains 
etc.). 



APPENDIX A. VARIATIONAL PRINCIPLE IN THE HEISENBERG MODEL 

In this Appendix we consider a variational approach to SSWT. This is a generalization of 
approaches of Refs. to quasi-2D case, which gives a possibility to improve the behavior 
of the magnetization at not too small J' / J. 

We apply the Feynman-Peierls-Bogoliubov variation principle for the free energy F = 
-TlnSp(e-^^) (see, e.g., |3|j5|) 

F<Fo + {H-Ho)o (115) 

where Hq is the trial Hamiltonian, Fq is the free energy corresponding to Hq, (...)o stands 
for the average with Hq. 

Consider first the case of a ferromagnet. Assuming the spin- wave character of spin 
dynamics we choose Hq as the Hamiltonian of non-interacting bosons and fermions 

Ho = —JolS + 7 E(^o - A)bib^ + {2S + 1)7 Jo E 
^ k k 

-/iEK^k+(25 + l)4ck] (116) 

k 

This expression differs from the Hamiltonian of the standard spin-wave theory by the factor 
7 which is a variational parameter describing the renormalization of the spin-wave spectrum, 
and by the presence of the Fermi operators taking into account the kinematical interaction 
between spin waves. Then we obtain 

^y^l^l^exp(-E^T2 

and 
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q pq q 

( i/o )0 = 7^0(5 - 5) - 7 E ^q^q - E [^q + (2'5 + l)iV/] . 



(118) 



where the magnetization reads 



5=(5 + -)coth^---EiV< 



Ef 1 



(119) 



and Ef = 'yJo — fi. The equation for 7 is determined from the condition dF/d^ = and has 
the form 



^ J ON dS- 



E^-qiVp^ + E^iV.- 



(120) 



In the antiferromagnetic case we use the trial Hamihonian 
Ho = -^JojS + ^7E ['^o(ak«k + blbk) - Jk(ak6^k + ^k^-k) 

k 

+ {2S + l)7^oE(4ck + 44) - E [4«k + bibi, + (2^+ l)(cU + 44)] (121) 

k k 

to obtain 

dS 



2^ "^q^ ^ "^0^ 



pq 



^7 q 



E 4-q-^p^:7^ + E 4-^q 



(122) 



where 



^q= V'(^^0-/i)2-(74) 



(123) 



and 



^=(5 + ^)coth|f-E 



(124) 



In the 2D case (J' = 0) we have 



''-'^^ 57 



(125) 



and a similar result with the replacement A'^q — >■ Lq. Then the equations ( |12(J| ) and 
reduce to 
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-'0 q 

7 - ^ = ^ E 



(FM), 
(AFM). 



(126) 



These equations coincide with those of Sect [III A[ At small J' / J the values of integrals over 
the Brillouin zone are determined by the region of small quasimomenta q^^Qy < J' / J. Then 
it is possible to put 



J, 



p-q 



~ Jo, -^q ^ Jo 



(127) 



in this case, which again leads to Eqs. (|126|) . If J'/ J is not small enough, the approximation 
( |127| ) becomes invalid and it is necessary to use the equation ( |12(]| ) or ( |122| ) up to the limit of 
the cubic crystal, J' = J, where equality ( |125| ) is again satisfied and passing to ( p.26| ) becomes 
justified. The same situation is realized for hypercubic lattices of any dimensionality in the 
nearest-neighbor approximation. 



APPENDIX B. THE SOLUTION OF INTEGRAL EQUATION FOR THE 
RENORMALIZED VERTEX OF SPIN- WAVE INTERACTION 



The equation for the vertex ( |100|) has a degenerate kernel. We search the solution in the 
form 



$(k, k - q; p - q, p) = J(Aq - Bp)k - 2 JA(p, q) 

Then we obtain after some algebraic manipulations 

2 



A 



+2A {jr;, - T, 

2 

5o/7 -R + M, 
A(p, q) = A[l-r,(Apq + 5g2 



l + 2q(p-q)(r,-M,/g^) 



B 



;i28) 



(129) 



where we have introduced the correction to the magnetization R = {Sq — S)/'~f and the 
non-uniform susceptibility x^o 
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R = 

V = ■ 

AqO 

and Tq, Mg are defined by 



T 



(fs 1 



IJI72 Js<q, (27r)2s2 + A 



T 



(J^yj (27r)2(s2 + A)[(s-q)2 + A]' 



(130) 



T 



\J\^^J (27r)2(s2 + A)[(s-q)2 + A] "^^^ 
T r (Ps SiSj 1 



IJI72 



(27r)2(s2 + A)[(s-q)2 + A] 2 



:i3i) 



Calculating the integrals in ( |131| ) and retaining only terms which are logarithmically diver- 
gent at A — > yields 

1, 



2 - -iJlr 



zz 
qO 



(132) 



where x^o given by (|103| ), and finally we get 



R 



T 



In 



27r|J|72 Ai/2 

Combining the above formulas we obtain the result (|101|) of the main text. 



(133) 
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CAPTIONS TO FIGURES 



1. Calculated temperature dependences of SRO parameters in S* = 1/2 quasi-2D ferro- 
(right-hand side) and antiferromagnets (left-hand side). The values of J' / J stand at 
the curves. Dots correspond to 2D case {J' / J = 0) without including the pseud- 
ofermion contribution, long-dashed lines and circles to calculations from Eqs. (p!3[), 
(p!4|), (|20|) and (|29|) , ( ^2]) respectively for J' / J = 0.3. Triangles mark the ordering 
temperatures. 

2. Temperature dependence of the (staggered) magnetization 5 for 5 = 1/2 quasi-2D 
ferro- (right-hand side) and antiferromagnets (left-hand side) with different values of 
J'/ J. Short-dashed lines show the results without inclusion of pseudofermions, long- 
dashed lines in the presence of pseudofermions, and solid lines correspond to the ap- 
proximation of effective SRO parameter of Sect. |111 B| . For J' /J = 1 the solid and 
long-dashed lines coincide exactly, and for J' / J = 0.01 long- and short-dashed lines 
coincide practically. 

3. Temperature dependence of the SRO parameters of quasi-2D ferromagnets for different 
spin values in the approximation of effective SRO parameter. 

4. Temperature dependence of the chemical potential of the boson-pseudofermion systems 
in quasi-2D ferromagnets in the approximation of effective SRO parameter. 

5. Temperature dependence of the magnetization 5 for a 5* = 1/2 quasi-2D ferromagnet 
with J'/ J = 0.1 in the external magnetic field. 

6. Diagrams corresponding to the spin- wave interaction contribution to different quan- 
tities (a) one-particle Green function of SSWT (b) effective vertex in RPA (c) non- 
uniform RPA susceptibility (d) correction to (staggered) magnetization 6S = S — S. 
Simple and bold lines denote the bare and renormalized one-particle Green functions, 
point stands for the bare vertex. 
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7. Temperature dependence of the relative sublattice magnetization a{T) of K2NiF4 in 
SWT, SSWT, RPA approaches and expansion in 0{N) model as compared with 
the experimental data (circles). The RPA' curve corresponds to inclusion of the func- 
tion given by ( |112| ). Short-dashed line is the extrapolation of the result of 



1/A^-expansion to critical region (see Ref. p6|). 
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